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Abstract — Most methodologies for modeling Transmission Line 
(TL) based channels define the input-output relationship in the 
frequency domain (FD) and handle the TL resorting to a two- 
port network (2PN) formalism. These techniques have not yet 
been formally mapped into a discrete-time (DT) block model, 
which is useful to simulate and estimate the channel response as 
well as to design optimal precoding strategies. TL methods also 
fall short when they are applied to Time Varying (TV) systems, 
such as the power line channel. The objective of this paper is to 
establish if and how one can introduce a DT block model for the 
Power Line Channel. We prove that it is possible to use Lifting 
and Trailing Zeros (L&TZ) techniques to derive a DT block 
model that maps the TL-based input-output description directly 
in the time domain (TD) block channel model. More specifically, 
we find an interesting relationship between the elements of an 
ABCD matrix, defined in the FD, and filtering kernels that allow 
an elegant representation of the channel in the TD. The same 
formalism is valid for both the Linear Time Invariant (LTI) 
and the Linear TV (LTV) cases, and bridges communications 
and signal processing methodologies with circuits and systems 
analysis tools. 

Index Terms — Power Line Communication, Phone lines, Coax- 
ial cables, Lifting-trailing-zeros, Transmission Lines, Time- 
varying block channel. 



I. Introduction 

TODAY many Home Networking (HN) technologies are 
available to the consumer. The general consensus is that 
HN is moving beyond simple data sharing among PCs and 
their peripherals and that multimedia applications will be at the 
heart of HN products. Although wireless Local Area Networks 
(LAN) based on the IEEE 802.11 standard are the most popu- 
lar HN solutions, they often suffer from poor RF propagation 
and from mutual interference |fl~). Hence, "wired" in-home 
(IH) connections are gaining new momentum, and high-speed 
HN technologies over IH power lines (PLs), phone lines (PHs), 
and coaxial cables (CX) are being intensely investigated. 
Among them, Power Line Communications (PLCs) represents 
the most accessible as well as the most challenging medium 
0, O. Initially, the interest for PLCs was targeted to using 
the existing outdoor PL infrastructure for delivering broadband 
Internet access [4|. More recently, PLC technology is emerging 
as an excellent candidate for many other applications, e.g. 
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for in-home, in-office and in-vehicle LANs, for smart grid 
applications, and in building control and automation networks 
0, 0, Q. These technologies are also being standardized 
in IEEE and ITU-T. PLCs are the scope of the recently 
confirmed IEEE P1901 baseline document [8], whereas ITU- 
T has recently consented to Recommendation G.9960 (also 
known as G.hn) for a unified IH networking technology 
operating over of all types of IH wiring |9]. 

TL-based channels like PLs, PHs and CXs have been 
traditionally modeled in the FD representing section of ca- 
bles and various common discontinuities like bridged taps, 
transformers, series impedances, splitters, etc., via Two-Port 
Networks (2PN) and their corresponding ABCD or Transmis- 
sion Matrices (TMs) $M, flU, EG), G3- While the modeling 
of CX and PH is today fairly well established, the modeling of 
PL is still under intense debate. The recent results published 
in lTT2l . fl3l . confirm that the ABCD-based formalism is 
adequate to describe signal propagation over PLs and that this 
conclusion holds also when grounding is present. However, a 
peculiar aspect of the PL channel that differentiates it from 
the other two IH TL-based channels is that the PL channel 
is a Linear and Periodically Time Varying (LPTV) channel as 
demonstrated by Canete et al. in lTT4l (see also Sect. [TT] for 
more details). 

Time-varying DT models for the TL-based channel are not 
available yet, since modeling has been traditionally carried 
out in the Laplace or Fourier Domain. These domains are 
inappropriate for time-varying systems. Some initial work on 
the modeling of LPTV nature of the PL channel can be 
found in |15|, fl6l . 03), ifTTl and references therein. These 
previous attempts at modeling do not provide a direct mapping 
between the modulated symbols and the output data, which 
typically are used to estimate as well as optimize the precoding 
method. Motivated by the fact that key advances in broadband 
wireless and DSL technologies were fostered by utilizing 
block transmission models and precoding strategies, our work 
contribution is twofold^ 1) we prove that block models similar 
to those used in wireless and wireline DSL channels can be 
used in the PLC conext as well; 2) we provide a DT model 
capable of handling the LTI, LTV and LPTV cases under the 
same formalism. The block transmission model we propose is 
immediately useful for the analysis of present designs, since 
the IH networking standards use schemes based on OFDM or 

1 Some initial results were presented at the 2008 IEEE ISPLC conference 
[18]; however, only the simpler case of Trailing Zeros was addressed there 
and the Zadeh decomposition was not given. 
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Wavelet-OFDM J8), ED, or single carrier schemes based on 
block DFE lfl9l 

The novelty of the method proposed here is that it addresses 
cases not contemplated by previous work on the modeling of 
the TV nature of the PL channel. For example, [ 14] allows the 
modeling of TV cases where the time variability is such that 
it can be modeled as a succession of LTI states thus allowing 
the use of a Fourier basis, whereas the goal of [17] is limited 
to the estimation of upper and lower bounds on the channel 
impulse response variations due to TV loads. Our work aims 
at mapping the modulated input signal to the output, modeling 
the variability of the channel using integral equations and then 
lifting and not a Fourier basis. As we will show here, this 
approach naturally leads to the design of precoders and bit 
loading methods, as well as facilitating receiver analysis and 
optimization. 

The paper is organized as follows. The transition from a FD- 
based to a TD-based model for TL-based channels is addressed 
in Sect. M for the LTV and LPTV cases. The lifted block 
representation of the DT model and the final input-output 
relationship of a TL-based channel is derived in Sect. [TTTJ in 
this section, we also verify the finite support of the employed 
LTV/LPTV/LTI filtering kernels and we propose an MMSE 
kernel estimator. In Sect. [IV] we show under what conditions 
the useful Chain Rule valid in the FD holds also in DT lifted 
form. Finally, we give conclusive remarks in Sect. [V] 

A. Notation and Acronyms 

In the rest of the paper we use the acronyms listed in 
Table [I] and use the following mathematical conventions: 
boldface upper (lower) case letters denote matrices (vectors); 
P x P channel matrices is denoted by the symbol H and 
its tall counterpart obtained by removing the last L columns 
is denoted by H; similarly, P x 1 vectors v[i] with trailing 
zeros are denoted as \r[i] = (vj [i], 0, . . . , 0) T ; element (fc, n) 
of matrix H is denoted as {H}fe.„; inverse and pseudoinverse 
matrices are denoted as _1 and respectively. 

II. From Frequency Domain to Time Domain: the 
LTV and LPTV Cases 

A general result of TL theory is that every uniform TL can 
be modeled as a Two-Port Network (2PN), thus allowing us 
to replace a distributed parameter circuit with a single lumped 
network. A 2PN has an associated 2-by-2 TM whose four 
elements (A, B, C, and D) are complex functions of the 
frequency. The steady-state relationship between current and 
voltage (in the FD) at the two ports of a 2PN is tied by T/, 
the forward TM (see Fig. [TJ: 
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In general, a TL-based link is made of several sections; each 
section consists of segments of different cables of various 
lengths spliced together. Series/shunt impedances, bridged 
taps, and other discontinuities may be present along the line. 
The TM of the overall link is simply obtained by multiplying 
the respective TMs of each portion of the network (Chain 



Rule). The TM formalism was first used in digital commu- 
nications applications by J.J. Werner for PH modeling |10|. 
Later it was adopted by Galli and Banwell for PL channels 
with and without grounding fl3l . 120 1. and by Chen for the 
CX case El. 
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Fig. 1: Schematic representation of a two-port network (2PN). 



It is useful to express output quantities as a function of input 
quantities inverting ([!}: 
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where the backward TM is T b = (T/)" 1 . 

The transfer function of the 2PN, i.e. the ratio of the voltage 
on the load Zq to the source voltage can be easily 
expressed in terms of ABCD parameters: 
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'0 TL,T,Ji *o ~ -^*o 
For some notable cases, the elements of the forward TM can 
be written in closed form: 

'Cable of length I: A = D = cosh( 7 (/)0; 

B = Z o (/) S inh( 7 (/)0; 

C = smh( 7 (/)Z)/W) 
A = D = h B = 0; C=l/Z 
A = D = 1; B = Z; C = 

(5) 

where -f(f) and Z Q (f) are respectively the propagation con- 
stant and the characteristic impedance of the cable. The behav- 
ior of 7(/) and Z a (f) for two typical IH wires is shown in Fig. 
|2j 7(/) grows very quickly with frequency whereas Z a (f) can 
be considered approximately constant in the High Frequency 
(HF) band and above. Furthermore, det(T/) = det(Th) = 1, 
i.e. the 2PN models of TLs and of shunt/series impedances 
are reciprocal 2PNs. In this case, ([T]i becomes: 



Shunt impedance Z: 
Series impedance Z: 



Vout(f) = D(f)V ln (f) - B(f)I ln (f) 
lout(f) = -C(f)V in (f) + A(f)I m (f) 



(6) 



Starting from the above FD representation of a TL, we will 
now introduce our DT models starting from the general LTV 
case and then particularizing it to the LPTV case. Then in Sect. 
Ill we will introduce the proposed block form representation 
(lifting) and show how the LTV, LPTV and LTI cases can be 
all handled under the same formalism. 
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TABLE I: List of acronyms used in the paper 



A r*vc\T\\Ttin 
1 U11V 111 


IVlCtllllllg 


l uiiy ill 


lVlCtllllllg 


2PN 


Two Port Network 


L&TZ 


Lifting and Trailing Zero 


AC 


Alternate Current 


LPTV 


T * 1 T~-\ l ■ -|-| m" T 7 - 

Linear and Periodically Time Variant 


CX 


Coaxial Cable 


LTI 


Linear Time Invariant 


DT 


Discrete Time 


LTV 


Linear Time Variant 


FD 


Frequency Domain 


OFDM 


Orthogonal Frequency Division Multiplexing 


FT 


Fourier Transform 


PH 


Phone Line 


FIR 


Finite Impulse Response 


PL 


Power Line 


HN 


Home Networking 


PLC 


Power Line Communication 


IBI 


Inter Block Interference 


RMS-DS 


Root-Mean-Square Delay Spread 


ICI 


Inter Carrier Interference 


TD 


Time Domain 


IFT 


Inverse Fourier Transform 


TL 


Transmission Line 


IH 


In-Home 


TM 


Transmission Matrix 


IIR 


Infinite Impulse Response 


TV 


Time Varying 


LAN 


Local Area Network 


TZ 


Trailing Zero 



(BLK): Bolden 124A Quad (solid): Coax RG-6 (dashed) 
Z (I) (RED): Belden 124A Quad (solid): Coax RG-6 (dashed) 




Fig. 2: Plot of the absolute value of f(f) and Z (f) versus frequency for 
two typical IH wires. 



A. The LTV Case 

In the LTV case |6]l does not hold because TV systems 
exhibit integral input-output equations in the frequency domain 
as well as in the time-domain. In fact, since the Fourier basis 
is no longer the basis of eigenfunctions in the case of LTV 
systems, we argue that there is no particular advantage in 
analyzing the system in frequency versus time. Furthermore, 
since communication signals are band-limited, a discrete time 
description is appropriate and sufficient to describe the input 
output relationships. Our objective is to lay down DT models 
that are flexible enough to handle the LTV, LPTV, and LTI 
cases. 

Since we are in the linear regime, in general we can write: 



Vout{t) = £v[Vin(t),kn(t)] 
iout(t) = £i[Vin(t),iin(t)] 



Hence, it is still possible to view the electrical quantities on 
Port 2 as filtered versions of the quantities on Port 1, where 
now the filters are represented by TV kernels, analogous to 
the time -varying impulse response introduced by Bello \22\, 
as follows: 



Vout(t) 



iout{t) 



Vin(T)d(t,t - T)d,T 

oo 

Un{T~)b(t, t - t)cLt 

-oo 
x> 

v in (T)c(t,t- r)dr 

oo 
oo 

+ I iin{T~)a(t,t - r)dr, 



(8) 



(9) 



(7) 



where d(t,t-r) = C v [S(t-T),0], -b(t,t-r) = C v [0,5(t- 
t)], -c(t,t - t) = Ci[S(t - t),0], and a(t,t - r) = 
£j[0, S(t — r)]. Our work assumes that an expression for the 
above TV kernels is available, but we point out that it is not 
trivial to obtain TV kernels in complete generality. Finding 
an appropriate mathematical representation is, however, a first 
essential ingredient to obtain the kernels either experimentally 
or through circuit analysis. 

There are differences between time variations caused by non 
uniform TLs and TV loads with uniform TLs. In |23|, |24|, the 
authors studied the case of non uniform TLs where the time 
variations are due to changes of the primary parameters of the 
TL. This case is seldom of interest since the most common 
TL based communications channels owe their time variability 
to loads only. The case of greatest interest in the context of 
data communications is when the time-variabilty of the loads 
are TV and the time-variability is periodic (LPTV). 

B. Zadeh 's Expansion Approach for the LPTV Case and Its 
Extensions 

The IH PL channel measurements in [ 14| give insight on 
how many significant harmonic components are present in the 
PL channel. lfT4l reported a median Dopppler spread of 100 
Hz, a 90%-percentile of 400 Hz, and Doppler components 
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up to B\ 



= 1, 750 Hz; thus, the coherence time of 



the channel is around I/-B5 



(max) 
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600/is. Interestingly, all 



observed Doppler components were quantized at multiples 
of the fundamental frequency of the mains AC cycle, i.e. 
fo = I /To — 50Hz. This confirms that the PL channel 
can be modeled as an LPTV channel. The LPTV nature of 
the PL channel has already been experimentally confirmed in 
02), lfl4"l and this paper is not concerned with confirming 
the validity of the LPTV model. Rather, our objective is 
that of finding a convenient representation and DT model for 
the overall channel response for the design and analysis of 
communication techniques over the PL channels. 

An LPTV system y(t) = C[x(t)] with period T is such that 
its time-varying impulse response £[5(t — r)] = h(t, t — t) is 
periodic with respect to time, i.e. h(t, £) = h(t + kTo, £), Vfc £ 
Z. In [25], Zadeh introduced the following expansion: 



h(t,r) 



E 

m= — 00 



(10) 



where the so-called harmonic impulse responses are defined 
as below: 

Mr) = ^ P h{t,T)e-^ mfot dt (11) 
J o Jo 

In Zadeh's expansion an LPTV channel is equivalently repre- 
sented as a bank of LTI channels whose outputs are modulated 
by Fourier harmonics with frequencies that are integer multi- 
ples of the fundamental frequency fo (see Fig. [5J: 
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\ ^ e j27rmf t 
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h m (r)x{t - T)dr (12) 
e i2 ™ /ot {/i m (r)«(t-T)} (13) 



The expansion given in ( fTT) uses the orthonormal Fourier 
basis {1/ y/%e j2mnfot / T rectTit)}™^^. For such an or- 
thonormal basis, if only M terms are used for estimating 
the channel response, then the channel estimate h(t, r) = 
Em=-M h m (T)e j2mnfot incurs a truncation error that 
is proportional to the energy of the truncated terms: 



E 

\m\>M 



\h m (T)\ 2 dr. 



(14) 



Remark 1: It is desirable to select a basis for which the 
summation in ( fT4] > grows slowly, i.e. a sparse basis. In 
the case of relatively abrupt discontinuities, the truncated 
terms of the Fourier basis in ( fTT| die out slowly. In these 
cases, one may replace Zadeh's model, with a more gen- 
eral basis expansion model for the TV response in the 
period ftr (t, t) — h(t,T)rectT (t), and obtaining back 
the h(t,r) = J2^L-oo hr {t - pTo,r). The expansion of 
h>T (t, t) is analogous to ( fT0| ) and obtained replacing with a 
Wavelet basis {p m {t)}m=-oc the orthonormal Fourier basis 
{l/V^' 2 ™ / o*/ T recfr r (i)}~ = _ 00 . The advantage is that, if 
appropriately chosen, the basis {p m (t)}^ = _ 00 can concentrate 
the energy of the channel response in fewer coefficients 




Fig. 3: Series expansion of an LPTV channel. 



h m {T) 126 1 . For an orthonormal expansion, the h rn (T) would 
be simply replaced by the following responses/coefficients 
h m {r) — J Q h(t, r)p^(t)(it. Seeking alternatives to the 
Zadeh's expansion is useful to reduce the truncation error that 
is inevitable in estimating or simulating the channel through 
these models. 

Remark 2: Note also that the presence of Dopplers in 
Zadeh's model for an LPTV system suggests that multicarrier 
transmission schemes will be inevitably affected by inter- 
carrier interference (ICI), regardless of the duration of the 
guard interval. ICI will cause greater performance degradation 
in conventional or windowed OFDM than in Wavelet-OFDM 
1271 or OFDM/OQAM |28| since Gaussian shaped filters in 
filterbanks are more effective at rejecting ICI. 

Several methods have been reported to model LPTV 2PNs 
(see HD, ED, ED). The first two generalize the ABCD TM 
formalism, considering LPTV responses, which have, however, 
no memory. The third reference ||3D derives the general 
response to an exponential waveform of a network that has 
periodically varying elements with the intent of extending the 
traditional spectral analysis methods to the LPTV case. More 
recently 11321 proposed to resort to a signals' sub-band decom- 



position, in the form x(t) = Y^t=-oo Xk(t)&' 2T 'T t ! along with 
Zadeh's filterbank expansion of the network response [25 1. 
Overall the authors express the linear relationship between 
the input and output subband frequency components using 
an infinite size Toeplitz matrix (called Harmonic Transfer 
Matrix) whose elements are i? ra _ TO (/ + m/o), where H m (f) 
is the FT of Zadeh's h m (t). Furthermore, [32| also presents 
an algorithm for the generation of symbolic expressions for 
the harmonic transfer functions of an LPTV systems starting 
from a system model in the form of a block diagram. What 
is inconvenient about this expansion, is that it requires the 
expansion of x(t) as ^fe^oo Xk(t)e : ' 27r T t and the Xk(t) have 
no particular relationship to the encoding and modulation of 
the signals. 

In general, a key advantage and difference of our represen- 
tation is that we combine different bases expansions for input 
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and channel. In fact, the basis chosen for the expansion of the 
LPTV channel does not in general have to be used as the basis 
of representation for the input signal as well, which is naturally 
represented through its samples, given that it is a band-limited 
communication signal. Furthermore, due to receiver filtering, 
the received data are also better represented through their 
samples, and many of the off bandwidth terms of the Zadeh's 
expansion are eliminated from the signal through filtering. 
Hence, the main contribution of this paper is to provide a 
new representation for the channel, that is consistent with the 
LPTV model but that, compared to the previous literature on 
the subject, operates directly in the sampled time domain. 
Our representation is centered on the idea that the memory 
of the LPTV PLC channels is approximately finite and that 
we can conclude that the harmonic responses are also of finite 
memory. Note that is is a sufficient but not necessary condition 
to have the overall LPTV PLC channel have finite memory. 

We verify this hypothesis in Section |III-B| relying on the 
fact that time-variations happen at a slow time scale com- 
pared to the ordinary symbol rates of PLC modems. That is 
reasonable, since the LPTV behavior is due to the fact that the 
electrical devices plugged in outlets (loads) contain non-linear 
elements such as diodes and transistors that, relative to the 
small and rapidly changing communication signals, appear as 
a resistance biased by the AC mains voltage. AC adaptors 
in particular, are highly likely to exhibit an LPTV input 
impedance. In fact, the periodically changing AC signal swings 
the devices over different regions of their non-linear I/V curve 
and this induces a periodically TV change of their resistance. 
The overall impedance appears as a shunt impedance across 
the "hot" and "return" wires and, since its time variability 
is due to the periodic AC mains waveform, it is naturally 
periodic. Furthermore, electrical devices are noise generators 
and, in view of Nyquist theorem, noise also appears to be 
cyclostationary. 

Although in all cases of practical interest in PLs one can 
exploit the above mentioned time scale difference, Zadeh's 
expansion for the entire channel response is useful in a more 
general set of cases. In fact, it allows for the development 
of adaptive techniques for the estimation of the harmonic 



responses of the channel (see Sect. III-Ai 



III. Lifted Representation of PL Channels 

Let us assume that a link is constituted of N cascaded 
sections and let hP-'"' ,N \t, r) be its TV overall channel 
impulse response. The signal at the output of the channel can 
be written as: 

/oo 
v a (T)h^-' N \t,T)dr. (15) 
-OO 

In DT, the equivalent relationship is: 

oo 

v ou t[k]= v.[n]h^-'^[k,k-n], (16) 

n— — oo 

where the expression of h^ 1 '"' ' N '[k,l] is given in ( fT7| ) at the 
top of the next page. p(t) is the transmit pulse shaping filter, 
andp'(i) is the receive filter and h^ 1 '"' ' N ^(t, r) is the baseband 



or baseband equivalent response of the LTV channel. Now, 
omitting superscript (1, • • • ,N) for brevity, our objective is 
to cast ( [16] ) in the following vector matrix model 11331 : 



w out [i] = Hi. v s [i\ + H M v s [i - 1] 



(18) 



This is possible by lifting the streams v ou t[&] and v s [n] 
into blocks as follows: {v s [z]}fc = v s [iP + k], {v out [i]}fe = 
v ut[iP + k], and by introducing the following two matrices 
(fc,ra = (),••■ ,P-1): 

{K i , }k,n = h[iP + k ) k-n} (19) 
{B i , 1 }k, n = h[iP + k,P + k-n]. (20) 

Matrices Hi,o and Hj i are banded and their structure is shown 
in Fig. |4] The underlying assumption for the representation in 
( |T~8] > is that the channel h[k, I] has finite memory L smaller 
that the block length P, i.e. P > L and h[k, I] = for I < 
and I > L. This allows us to limit IBI only to two consecutive 
terms. If this condition were not met, then the output block 
v ou t[«] would be a function of infinite IBI terms. 

Since only the top right corner of H,;.i is non zero, having L 
zeros in the last L entries of every size-P transmitted block al- 
lows us to eliminate completely IBI thus yielding to simplified 
expressions (TZ case). In fact, posing v s [i] = (v s [z], 0, . . . , 0) 
allows us to replace ( fl8] > with an IBI free counterpart: 



v otlt m 



H 



i,ov s [i\ 



HiV s [ 



(21) 



L+1 



H l0 = 





Hu = 



Fig. 4: Time-varying block channel matrices have a lower-banded (left) and an 
upper triangular (right) structure of order L. The elements of these matrices 
are given in (19) . 

Remark 3: Even though H^rj in ( f2T| i is a lower triangular 
matrix and, therefore, it is invertible, its condition number is 
significantly affected by the channel response. Such a matrix 
is well conditioned only when the channel is minimum phase. 
Interestingly, its tall counterpart Hi, obtained by removing 
the last L columns of H^, is typically well conditioned and 
its properties are discussed in detail in (34]. In general, in 
solving for the input, one should invert the over-determined 
system JTx[i] rather than inverting Ho, to avoid numerical 
instability. 

A. Estimation of Zadeh 's Harmonic Impulse Responses 

Let us assume that /o <C W, where W is the 
bandwidth of the receive filter p'(r), and let h m [l] = 
S-oo I-oo h m (T)p(lT s - t - 0P'(0 dTd $< where h m(T) are 
the baseband or baseband equivalent harmonic responses in 
( [TTj ). Note that, given the receiver filtering, any component of 
the LPTV channel outside the receive filter will be filtered out. 



S. GALLI et ah: DISCRETE-TIME BLOCK MODELS FOR TRANSMISSION . 



oo /»oo 



hP'~> N )[k,l]= / / h^-^(kT s -^r)p(lT s -r-0 P '(Odrd^ (17) 



CX3 •/ — OO 



Introducing the following matrices for k,n = 0, • • • , P — 1, 
and I = 0, 1: 

{H m ,/} fc! „ = /i m [ZP + fc-n] (22) 
{n m } fc ,„ 4 e ^™/° feT .,5[n-fc], (23) 

Zadeh expansion for the channel matrix in block form is as 
follows: 

{H,.;k, - ^ ft m [/P + A - n ]e^ ro '°( iP + fe ^ (24) 

= ]T e j2 ™ /o4PTs {fi m H m ,i} fc ,„. (25) 

?n 

In the case of TZ, assuming that the h m [I] have finite memory 
L, denoting by h m the vectors of their non zero coefficients, 

and by {Hi} k , n = T, m h m[k - n }e^^ m( - w + k ^ for k e 
[0, P — 1], n G [0, P — Z — 1] we can analogously remove the 
IBI and have: 

Vout[i\ = ^2e^ m ^ iPT 'Sl m H m y s [i) (26) 

m 

= ]>> j2 ™ /o * PT ^ m <I>(vJ i ])h m , (27) 

m 

where the P x L Toeplitz matrix <&(v s [i]) has first column 
v s [i] and first row (v s (iP), 0, . . . , 0). 

Assuming that the response can be well 
approximated with a limited number of harmonics 
m G [— M, M], stacking all the kernels in 
h = (h^ M , . . . , hjy-) T and forming the matrix *(v s [i]) = 
(e-^ M /° iPT *n_ M *(v s [i]), . . . , e^ M /° iPT *fi M *(v s [i])), 
this means that: 



v out [«] = *(vji])h 



(28) 



and a unique MMSE estimate of Zadeh's harmonic impulse 
responses exists if the matrix is at least square (P > (2M + 
1)L) and invertible. The MMSE estimate of h is then: 

h = * t (v s [ l ])v OMt H. (29) 

Note that estimate ( f29) is more accurate than the one proposed 
in equation (4) of |14| which is an ML estimate that neglects 
noise and also requires the use of a Fourier basis. 

We have just shown that ( |2"9| ) allows to estimate the 
LPTV channel from the sampled observations at the receiver. 
The error in the estimated channel response can be readily 
quantified after recognizing that it includes three terms: 1) 
the residual MSE due to noise in the observations which is 
given by trace((^ H (v Ji])R n \J r (v s [i])) _1 ), where R„ is the 
noise covariance matri:?(rl 2) the series truncation errors given 
in ( |T4] > due to considering a limited number of harmonics 
hi € [— M, M\ and, last but not least, 3) the approximation 
of the responses with Finite Impulse Response (FIR) discrete 



2 The noise covariance matrix R,j is seldom diagonal as noise in PLs is 
colored (see |2| and references therein). 



time filters. The nature of this last approximation will be 
discussed and quantified in the next Section. 



B. On the Finite Support of the ABCD Kernels 

The time scale of the excitation given by the AC signal is 
very slow compared to both the duration of the responses of 
the LPTV 2PN as well as the Nyquist interval for transmission 
over PLs, given the broad bandwidth of the input signal itself. 
This, we argue, implies that not only the load equations can 
be linearized, but also that their response can be approximated 
resorting to the solution of the small signal equivalent circuit, 
as if it were stationary. Hence, it is possible to map the ABCD 
parameters in linear impulse responses. Based on |5]l it is 
clear that the kernels a(t, r), d(t, r), which are different from 
one only for an LTI section of cable, are always LTI impulse 



responses a(r),d(r) (see Section III-B1 1 while, for series and 
shunt impedances, this yields: 



b(t,r) 



c(t,r) 



+ 00 

E 



b m {r)e 



j2Trmf t 



+oo 

E ■ 

in— — oo 



i( r )e 



j2irmf t 



(30) 



(31) 



The above expressions are TV once the relationship between 
the periodic AC bias and the small signal circuit parameters 
that determine their impedance is made explicit. Given the 
Z(t, f) the Zadeh expansion terms are as follows: 

f+oo r T 

Z(t, f)e-^ m}at e 32lTST dtdf (32) 



,(t) 



1 

% 
1 



-OO 



To 



-j2Trmf t j2rrfT 



dtdf (33) 



Note that the low pass equivalent kernel, which are of greater 
interest for PLC, can be obtained by simply replacing Z(t, /) 
with Z(t, f + / c ), i.e. down-shifting the expression of the 
impedance by the carrier frequency f c . 

At this point, the question we want to address is whether 
an analytical path towards generating the Hjj on the basis of 
the TL formalism exists or not. Specifically, the issue at hand 
is whether or not one can resort to a DT lifted representation 
of the TV kernels in the TL equations (jTOf to get to ( fT8j ). 



As mentioned in Sect. Ill the underlying assumption for 
exploiting lifting is that the kernels must have finite memory, 
and the objective of this section is to verify it. 

Based on the separation of time scales discussed earlier, in 
our analysis we treat each section of the PL link as LTI. When 
representing LTV loads, this is exact for the parameters A and 
D and it is valid for parameters B and C when the separation 
of time scales holds. When representing sections of cable or 
bridged taps, this will also be exact as they obviously represent 
LTI channels (assumption of uniform TLs). 



6 



SUBMITTED TO THE IEEE TRANSACTIONS ON COMMUNICATIONS, JUNE 2009. REVISED NOV. 2009 



For the LTI case, one can rewrite |6]) expressing the electri- 
cal quantities on Port 2 (output) as the result of a convolution 
between the electrical quantities on Port 1 (input) and LTI 
filtering kernels: 



Vout(t) = 


/ v in (r)d(t- 

J —OO 


r)dr 


(34) 




poo 

- / iin(r)b{t 

J —oo 


— r)dr 




iout{t) = 


poo 

- / v in (T)c(t 

J — oo 


— r)dr 


(35) 




poo 

+ / Hn(T)a(t 

J —oo 


— r)dr 


(36) 



where a(t), b(t), c(t), and d(t) are the impulse responses of the 
2PN filtering kernels and they are obtained by Inverse FT (IFT) 
of A(f), B(f), C(f), and D(f), respectively. Relationship 
( |36"1 > holds for both the cases where the filtering kernels 
represent a single 2PN of the kind of <|3j or the overall TM 
of a cascade of 2PNs. 

To verify that a(t), b(t), c(t), and d(t) have finite memory, 
we study the following three cases: a single section of cable, 
single and cascaded shunt impedances, and single and cas- 
caded series impedances. Without loss of generality, we resort 
to the baseband kernel model. 

1 ) Single Section of Cable: The forward ABCD parameters 
have the closed form expression in (151. A single section 
of cable acts on the transmitted signal as a low-pass filter 
with decreasing cut-off frequency as the length of the section 
increases; thus, it can be modeled as an IIR-like filter and, the 
longer the cable, the longer the impulse response and the more 
the IIR-like behavior will be. This can be verified analytically 
by noting that filters l/A(f) and 1/D(f) have the following 
known FT pair [35 1: 

_2 



sech(at) — sech 
a 



(37) 



and limt_ i . 00 sech(at) = 0. It is immediate to recognize that 
the DT version of the inverse filters of A(f) and D(f) are 
likely to be well approximated by FIR responses while the 
impulse responses a(t) and d(t) are better approximated by 
IIR filters and thus would result in IBI from infinite blocks. 
However, truncating the IIR response leads to a small error 
when the section of cable is not too long. This can be verified 
looking at Fig. [5] where we can see that the filtering kernel 
a(t) =IFT{A(/)} for a single section of telephone AWG 24 
cable is FIR-like until lengths of 1.5 kft but starts becoming 
IIR-like around 2 kft. For example, at a length of 1.5 kft 
we have ascertained that the energy of a(t) truncated after 
0.75/js is equal to around 98.5% of the energy of the non- 
truncated response. We also note that a(t) is equal to the 
impulse response of the link for the lengths shown in Fig. 
|5](a). This fact is not surprising as, when the second port is 
unterminated, the output of the filter A(f) is equal to the 
voltage on the second port and the ratio of the voltages at the 
two ports is equal to the impulse response for an ideal source. 

For the case of outdoor topologies, one can split the longer 
cable sections in multiple sub-section of smaller lengths, find 
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Fig. 5: Plot of the kernel a(i) =WT{A(f)} for the case of single section of 
an AWG24 cable of length I feet, on a band between DC and 30 MHz after 
raised cosine filtering with roll off equal to 0.5. (a) I = 100, 500, 1000, and 
1, 500, from left to right; (b) I = 2, 000. 



the DT lifted version of these sub-sections and then consider 
the cascade of these multiple sections to model the whole 
cable. In the FD, this is easily done by using the Chain 



Rule. In Sect. IV we will show that it is always possible to 
handle tandem connection of 2PNs even in the DT lifted form; 
furthermore, we show how a DT lifted equivalent of the FD 
Chair Rule exists under the TZs assumption. 

As an alternative for the case of cable sections of longer 
length, we can also introduce new kernels rewriting ([TJ as 
follows (A(f) and D(f) are always non-zero): 

B(f) 



V ou t(f) 
lout(f) 



1 
1 



V in (f) 



A(f) 

-c(f) 



lout(f) 

V out (f) 



(38) 



D(f) D(f) 

The new filters a(t), (3(t), 7(i), and £(£) have impulse 
responses that are obtained by IFT of 1/A(f), -B(f)/A(f), 
1/D(f), and —C(f)/D(f), respectively. The transfer func- 
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tions of these new kernels are of the kind sech(/) and tanh(/) 
Since the following FT relationships hold 

J sech(ai) 
1 cosech(ai) 



I sech(^/) 
4»£tanh(£/) 



(39) 



and we also have that lim t _ i . 00 cosech(ai) = 0, these kernels 
decay exponentially. The filtering kernels a(t), f3(t), j(t), and 
£(t) for the case of a 2 kft AWG 24 cable are shown in Fig. 
|6](a), and a(t) is compared to the impulse response h(t) of the 
link in Fig. |6](b). Besides verifying that the filtering kernels 
are well approximated by FIR filers, we also note that a(i) 
coincides with h(t) for the first (apart from a scaling 
factor) and then is followed by echoes of alternating sign. 



0.02 



x 10 

20 r 



E 5 - 



40 60 
Time (us) 



(a) 



12 15 18 21 24 27 30 33 36 39 
Time (us) 



6p,s from each other which corresponds to a signal traveling 
4,000 feet, the roundtrip of the length of the cable. This 
result is in agreement with ll36l where a model for the echoes 
bouncing between discontinuities was developed for the PH 
case. Although infinite echoes are theoretically present, only 
few echoes are non-negligible due to cable attenuation. The 
fact that IFT{A(/)} and IFT{1/A(/)} can well approximate 
the impulse response of a TL link has never been noted and 
it is reported here for the first time. 

2) Single and Cascaded Shunt Impedances: Shunt 
impedances model appliances plugged into an outlet as well 
as bridged taps [ 10 1 . In the vast majority of cases, a bridged 
tap is composed of a single section of cable terminated on 
a constant impedance Z t . For shunt impedances, the kernels 
are a(t) = d(t) = 6{t), b(t) = 0, and c(t) =IFT{1/Z m (/)} 
(see (|5]l) where: 

A(f) ■ Z t (f) + B(f) 



Z in {f) = 



(40) 



C(f) ■ Z f (f) + D(f) 

The above expression can be simplified further in the HF 
region as follows (see Fig. [2j: 



c(f) 



(41) 



(b) 



where we have exploited lim^^oo cosh(a;) « linx^oo sinh(x) 
and posed lim/_ i . 00 Z a (f) — Z%°. An IH bridged tap is 
often terminated with an open, a short, a resistor, or a small 
capacitor, so that in the HF band and above C(/) is constant 
and we can write c(t) « kd(t) and its kernels have no memory 
at all. The same considerations apply to N cascaded shunt 
impedances since their combined TM has the same form of 
((5) but with C = J2? =1 l/Z®, where Z® is the i-th shunt 
impedance in the cascade. 

3) Single and Cascaded Series Impedances: We have that 
o(t) = d(t) = S(t), c(t) = 0, and b(t) is equal to the IFT 
of the series impedance (see |5]l). As for the shunt impedance 
case, it is easy to verify that the TM of the cascade of N 
series impedances still has the same form of |5]l but with 
B = z(i) > where z{i) is the *" th 

series impedance in 

the cascade. We assume here that b(t) is always FIR-like in 
all cases of practical interest. 

Remark 4: The considerations made here for the LTI case 
are novel, and the existence of a (approximately) finite support 
for the LTI kernels is reported here for the first time. Since 
LTV TLs are composed by LTI sections (wires) and LTV 
sections (shunt/series impedances), the capability of handling 
both cases under the same formalism is of interest. 



Fig. 6: Case of a 2,000 ft long AWG 24 cable, (a) Plot of the filtering kernels 
a(t) = 7(t), /9(t)/50, and 150CW of |HJ- ( b ) Comparison between the 
impulse response h(t) (black) and the kernel 4ct(t) (red). 

Again this is not surprising as the output of the filter 1/A(f) 
when the second port is unterminated is equal to the voltage 
on the first port and, therefore, these alternating echoes are 
the successive reflections between the unterminated second 
port and the first port. In fact, echoes are spaced around 



C. The TL Input-Output Relationship in Lifted Form 

To express the TL input-output as in ( [IS) , we will first tackle 
the problem of expressing in DT lifted form the continuous 
TD input-output relationships in ( fTO) . This, in turn, will give 
us a DT equivalent of the Chain Rule to handle the cascade 
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of 2PNs. The DT equivalent of ( flu) can be written as follows: 

oo 

v<wt[k]= E v in [n]d[k,k-n] (42) 

n— — oo 

OO 

ii n [n]b[k,k - n] 



[k] 



71— — OO 
OO 

n— — oo 

oo 



E 



\k,k- 



]a[k, k 



(43) 



where the LTV DT kernels above are defined as in (17) . Other 
relationships that will be useful for finding the DT input-output 
relationship of a TL link are based on Kirchoff laws: 

v in [k] = v s [k] ~z { Q s) i m [k] (44) 
v ut[k] = z^ L) i out [k] (45) 

Since we have ascertained the existence of either exact or 
approximated FIR filtering kernels that have finite memory L, 
we can rewrite the previous DT equations in lifted form for 

P > L: 



Vout[i] = T>i : Qv in [i\ + T) it iv in [i - 1] 
- B i)0 ii„[i] - B iA Un[i - 1] 

wH = -Ci,oVi n [i] - Ci t iVi n [i - 1] 
+ A ij0 im[i] + A M i in [i - 1] 



Voutli] 



(«)• 



(46) 
(47) 

(48) 
(49) 



where we have posed: {v s [z]} fe = t; s [iP + k], {v out [i]}k = 
v out [iP+k], {v in [i]} k = Vin[iP+k], {w[i]}k = w[iP+fc] 
and {i ln [i]} fc = i m [iP + fc]. 

The P-by-P matrices Aj ; *, Bi Cj * and Dj * are defined 
as in ([19]), e.g. (fc, n = 0, • •' • , P - 1): 

{Ai ti -j} k>n = a[iP + k,(i- j)P + k - n] (50) 

On the basis of the considerations made at the end of Sect. [H] 
we can give explicit expressions to some of the above matrices 
as shown in Table [TTJ 

Assuming the system is at rest before the first block of the 
excitation v s [1] is transmitted, one can use the four equations 
(|47)-(|49) to solve for the four unknown blocks v out [l], i out [l], 
v m [l], and i in [l] by posing v in [0] = and i in [0] = 0. At the 
next transmitted i-th blocks (i > 1), one must also carry the 
IBI terms as known vectors in the same four equations and 
solve again for the same four unknown blocks. At every step, 
a system of AP unknowns and 4P equation must be solved. 

If TZs are applied to the transmit signal v s [i], one can find 
the solution of the system without IBI. In fact, after eliminating 
Vi n [i] using ( |48| ), we can rewrite (|47]i-(|49) as follows: 



v out [i] 




D il0 


UmtW 







-(4 s) D ii0 +Bi, ) 



( z o Q 



(51) 



Now, by exploiting also d49| we can write: 



[ I 



I ] 



D 



i,0 



-(4 S) Dz,o- 



B 



8,0 J 



C;,o ( 2 o ^i,o + A.,-,o) 



(52) 

and solving for we obtain ( |53] l the following expres- 

sionL: 



where matrix S 



= 4 (D. ii0 + 4 L) Q,o)v s w 

is defined as 



(53) 



"i,0 — V Z U i,0 I -°i,0 T ^0 ^0 °i,0 ~r A i,oJ 

(54) 

In the above equation we have taken into account the numer- 
ical instabilities that may arise when inverting lower banded 
matrices as mentioned in Remark|3] Thus, rather than inverting 
matrices Aj.o, B^o, C^o, and D^o we have resorted to t heir 
tall counterparts A^ , Bi,o, C^o, and -D^o, respectively. We 
also point out that S^q resembles very closely the denominator 
of equation |3]l which yields the channel transfer function in 
the FD. 

Finally, substituting the expression for ij n [i] in ( f5T| l, we 
obtain the input-output expression given in ( f55| ) shown at the 
top of the next page, where in the last row we have replaced 
the Px 1 input signal vector v s [i] with TZs with the (P—L) x 1 
vector v s [i] defined as v[i] = (vj[i], 0, . . . , 0) T . 



IV. DT Lifted Chain Rule 

It is of particular interest to check whether it is possible to 
extend the ABCD Chain Rule defined in the phasor domain 
to the DT lifted case. Let us start by expressing the input 
output-relationship of a 2PN in the following matrix form: 



Vout [i] 






B, 








A, 



UnH] 



(56) 



where 



Aj = 


[A lj0 :A 4 ,i] 








Bj = 


[Bi^rB^i] 


and < 


Vi„[i] = 


_ v in [« - 1 


c i = 






Un[i] = 


lin [^] 




[Di, :Di,i] 






_ Un [i ~ 1] 



(57) 

We can then state the following Theorem. 

Theorem 1: The input-output relationship of the cascade 
ofN 2PNs in DT lifted form is: 

r(l,-,JV) 



v ou4 [i] 









C 



(1,-,7V) 



-bI 1 --^ 



Un [i] 



where 



.*) _ r A d,- 



B 



(1,-,JV) 



\B 



a,- 

8,0 



(58) 

(59) 
(60) 
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Vout[i] = [D it0 - (z ( s) Di, + B ii0 )sT (£> i)0 + z { L) C i)0 )]Y a [i] (55) 



(i,- : A'): r ,(V-,jv)j 



[qr"":c 



{^1,0 J' 



(61) 
(62) 



awaf w/zere f/ie sub-blocks of these matrices can be computed 
recursively as follows: 



C 
D 



(i,-,fe) 

»,0 



r ,(fe) T ,(l,-,fe-l) 
r .(fe)„(l,-,fc-l) 



AWcgj-'*- 1 ) (67) 



B^C^— (68) 



,fe-i) 



Proof: The proof follows easily from from Theorem [T] 

□ 



A 


• ,fe) 


r,(t) n (l,- ,fe-l) 
= Ho^i.O 


L A ( fe ) A C 1 .- ,k- 


-1) 




A (v 

1,1 


-,fc) 


_ n (fe) T >(l,-,'=-l) 




-1) 








' •**■»,! "i-1,0 


1) +A (fe) A (l,- 


,k- 


1) 




-,fe) 


_ n (fc)p»(l.-,fc-l) 




-1) 




B&' 


-,fc) 


- r»( fc )pi( 1 '- ,fc - 1 5 


_i_ r»( fe )'R( 1 '"' > fe 


-1) 








+ D !,l A i-1,0 


i) , B W A (V~ 


,k- 


1) 


c (i,- 




_ p(fc) n (l.- ,fc-l) 


, A« r (V" ,k 


-1) 




r (V 


-,fe) 




i ri(*)r»(i.— : fe 

+ U i,0 U i,l 


-1) 








, A (fc) r (l,-,*- 

"T" °i-l,0 


i) , »Wp(V" 




1) 


D& 


• ,fe) 


- nWnt 1 '™.*- 1 ) 


, T>(fc) r i(l,- ,fc 


-1) 






-,fc) 


- nMn' 1 ''"'*- 1 ) 


_i_ n(*0-n(V-- .ft 


-1) 








, r>(fc)/-«(l,— 




,k- 


1) 



Proof: Using d58l for N = k, the claim of the Theorem 
is obtained by recognizing that ' 

.(l,-,fc-l) u _ .(fc)r, 



v,-„ \i] and 



that i\ 



and that the product of two upper 



(AO 



triangular matrices is always zero, i.e. HQ x 



H^- 1 ' = 0. 



□ 



Remark 5: The chain rule valid in the phasor domain allows 
us to multiply directly the transmission matrices to obtain the 
transmission matrix of the overall link. We recognize that this 
is not possible in the lifted form when IBI is present. However, 
if we add trailing zeros to the input signal v s [i], then we 
obtain the DT lifted equivalent of the Chain Rule of the phasor 
domain as shown in the following Corollary. 

Corollary 1: If the input signal has L trailing zeros, then 
the Chain Rule in the DT lifted form has the following 
expression for the cascade of N systems: 



Vout [i] 









D <1," 



J i,0 



n 

k=N 



r (fc) 



_ B (1,» 

D !,0 
A (l,- 

R (fe) 
A (k) 



,N) 



1 i I 



1 ; i 



(63) 



(64) 



and the sub-blocks of matrix ( |63| > can be recursively expressed 
as follows: 

A (i,-,fc) _ r ,(fe) T ,(i,...,fe-i) , A (fc) A (i,-,fc-i) 



B 



H,0 

i,0 ~~ 



D (fe) B (l,- 



,k-l) 



. A W A (L- 

. B W A (1.- 



,fc-i: 



(65) 
(66) 



V. Conclusions 

We have proposed a DT equivalent model for TL-based 
communications channels which are usually modeled exclu- 
sively in the FD and only allow the treatment of the LTI case. 
The formalism adopted for our models allows us to address 
both the LTI and the LTV/LPTV cases which are of practical 
interest as they represent common home networking environ- 
ments. Not only our method is suitable for more realistic 
simulations when the TL is LTV/LPTV, but also naturally leads 
to design of precoders and bit loading methods that are tailored 
to better cope with the periodically TV channel distortion since 
they directly relate the modulated symbols with the received 
data. These designs are going to be the subject of future work. 
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